Elementary excitations in trapped Bose gases beyond mean field approximation 
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Using hydrodynamic theory of superfluids and the Lee- 
Huang- Yang equation of state for interacting Bose gases we 
derive the first correction to the collective frequencies of a 
trapped gas, due to effects beyond mean field approximation. 
The corresponding frequency shift, which is calculated at zero 
temperature and for large N, is compared with other correc- 
tions due to finite size, non-linearity and finite temperature. 
We show that for reasonable choices of the relevant parame- 
ters of the system, the non-mean field correction is the leading 
contribution and amounts to about 1%. The role of the de- 
formation of the trap is also discussed. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80. Pj, 67.90. +z 



The experimental measurements (T]-|j| of the collective 
oscillations of Bose condensed gases confined in magnetic 
traps have provided an excellent confirmation of the pre- 
dictions of mean field theory (see Q for a recent theoret- 
ical review). The accuracy of the mean field predictions 
is not surprising since in these gases the average distance 
between particles is significantly larger than the range of 
interatomic forces. Typically, the gas parameter n(0)a 3 , 
where n(0) is the density evaluated in the center of the 
trap, and a is the s-wave scattering length, is smaller 
than 10~ 4 . According to the theory of Lee, Huang and 
Yang (LHY) the first corrections to the mean field 



and to be 



predictions are expected to behave like V 
consequently of the order of 1% or less in these systems. 
While such corrections are too small to be observed in 
the density profiles or in the release energy, they might be 
observable in the frequency of the collective excitations 
where the accuracy of measurements is much higher. For 
example an accuracy of 0.3 ~ 0.4% has been already 
achieved in the experiment of || . 

Measuring effects beyond mean field theory is a chal- 
lenging task and would open new perspectives in the 
many-body investigation of these novel systems. So far 
the theoretical investigation of these effects has been lim- 
ited to the equilibrium properties, either including first 
quantum corrections in analytic form J^Q, or through 
numerical simulations based on Monte Carlo methods || . 
The purpose of this work is to provide an analytic calcu- 
lation of the first corrections to the collective frequencies 
due to non-mean field effects. These corrections are cal- 
culated in the large N limit and at zero temperature. 

Our starting point are the hydrodynamic equations of 



superfluids (see for example M) 
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holding at zero temperature. Here n is the density of 
the system, v is the velocity field and /i is the chem- 
ical potential. Equations (|l|||) permit to describe the 
low frequency collective excitations also in non-uniform 
Bose superfluids, provided the density profile varies on 
a macroscopic scale and one can use the local density 
approximation for the chemical potential 



fx(r,t) = fj,i(n(r,t)) + V ext (r) 



(3) 



where /J.i(n) is the chemical potential calculated for a 
uniform gas at density n and V ext is the external confining 
potential. 

In the following we will consider the linearized regime 
of eqs.([IJ) and (|[). We write n(r, t) = n(r) + Sn(r,t) and 
/i(r,t) = /io + 8fii(r,t) with Sui = (d[j,i/dn)5n, so that 
(Q) and (||) can be rewritten in the useful form 
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The ground state density n(r) entering equation ([|) can 
be easily calculated by imposing the equilibrium condi- 
tion fiQ = Hi (n(r)) + V ext (r) where \i a is the ground state 
value of the chemical potential, fixed to ensure the proper 
normalization of n(r). 

Equations (|l]-||) do not necessarily require that the 
trapped gas is weakly interacting. It is also worth notic- 
ing that the density n entering these equations should not 
be confused with the density of the condensate, which 
in general does not obey equations of macroscopic type. 
Only in the weakly interacting limit, where quantum de- 
pletion effects are negligible, can the density of the sys- 
tem be identified with the condensate density. In this 
case eqs. (|l[^) are equivalent to the time dependent 
Gross-Pitaevskii equations for the order parameter, eval- 
uated in the large N limit (see, for example, [|l0| ) . 

According to LHY theory, the chemical potential of a 
uniform interacting Bose gas is determined by the low 
density expansion: 
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where g = Airh^a/m is the interaction coupling constant 
and a is the s-wave scattering length. Eq.(^) represents 
a major result of many body theory and accounts for non 
trivial renormalization effects of the interaction coupling 
constant. It provides the first correction to the result 
fj, = gn given by lowest order theory, hereafter called Bo- 
goliubov or mean field approximation. The LHY equa- 
tion of state (||) can be derived starting from Bogoliubov 
theory. In this scheme the energy of the system, includ- 
ing the zero-point motion of elementary excitations, is 
given by E= \Ngn + |Y^ P7t0 [e(p) - V 2 /2m- gn] where 
e(p) is the energy of elementary excitations in Bogoli- 
ubov theory. The zero point energy contains an ultravi- 
olet divergency at large p which is cured by the proper 
renormalization g — ► g(l + g^ J2 P ^o m /p 2 ) of the cou- 
pling constant, so that one finally obtains a convergent 
result for the energy of the system, yielding result (|5|) 
for the chemical potential. It is worth noticing that the 
LHY correction does not involve additional parameters 
with respect to Bogoliubov theory, being fixed by the 
scattering length and by the density. 

Using the LHY equation of state (§), the equation for 
the ground state density can be solved by iteration and 
one finds the result [pi 



n(r) = n TF (r) - tw%(r) 



(6) 



where nrf(r) = (p.Q — V ex t(r))/g is the so called Thomas- 
Fermi result for the ground state density [p"lf and a — 
(32/3V^)a 3/2 - 

Notice that in the same scheme the condensate density 
n c is given by the expression n c (r) = riTF(r) — j&n?np(r) 
which is smaller than the density n(r) and yields H the 
result 



N 



vaMo) 



(7) 



for the quantum depletion of the condensate of an atomic 
gas confined in a harmonic trap. 

Using (Jsj) one can write (dfii/dn) = g(\ + 3/2anl/p) 
so that eq.(0) takes the form 



mui 2 5n + V (gi%T 'Sn) = ~l^ 2 (oigri^pSn 



(8) 



Equation provides the appropriate generalization of 
the zero-th order hydrodynamic equation 



mu) 2 5n + V {gn^F^Sn) = 



(9) 



used in |T(| to evaluate the collective frequencies in the 
large TV, Thomas- Fermi approximation. Eq.^) admits 
analytical solutions if the external potential is harmonic. 
In particular, for an isotropic trap (V ex t = \muj 2 r 2 ) the 
solutions of (||) obey the dispersion relation ]10[ ] 



uj(n r , i) = u (2n 2 r + 2n r l + 3n r + £) 1/2 . 



(10) 



where n r is the number of radial nodes and £ is the an- 
gular momentum of the excitation. Equation ( [l0|) shows 
that the collective frequencies in the mean field approxi- 
mation are fixed, apart from geometrical factors, only by 
the oscillator frequency, a quantity measured with very 
high precision in experiments. This is a remarkable fea- 
ture exhibited by these harmonically trapped gases which 
allows for a safe investigations of small corrections. 

Once the solutions of (||) are known, eq. (||) can be eas- 
ily solved by treating its right hand side as a small per- 
turbation. One finds that the corresponding frequency 
shifts obey the general equation 
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where Sn are the solutions of (g) and u> are the corre- 
sponding frequencies. The integrals of ( |TT| ) extend to the 
region where the Thomas-Fermi density is positive. 
In the absence of trapping the gas is uniform and 



the solutions of eqs. 



have the form Sr 



and exhibit a phonon dispersion uj = cq. In this case 
Eq.(||) (or, equivalently, ([n])) gives the Beliaev result 
jUI Sc/c = 8y/a 3 n/Tr for the shift of the sound velocity 
with respect to the Bogoliubov value c = yjgn/m, calcu- 
lated at the density n. Notice that even in the uniform 
case n differs from utf because of eq.(||). The shift of 
the sound velocity is consistent with the change in the 
compressibility mc 2 — ndu/dn associated with the LHY 
correction in the equation of state (^|). 

Eq.([il"l) shows that the so called "surface" oscillations 
5n = r e Yi m , satisfying the condition V 2 <5ri = 0, are not 
affected by the LHY correction. Physically this behavior 
follows from the fact that in the long wavelength limit 
the hydrodynamic surface oscillations are entirely deter- 
mined by the external field. For spherical trapping these 
solutions obey the dispersion law u> = y/tuiQ which is 
simply obtained setting n r — in (|To|) . 

In order to observe effects beyond mean field one has 
consequently to focus on compressional modes, which are 
sensitive to the equation of state. The lowest mode in 
a spherical trap is the monopole (breathing) oscillation 
(ri r = 1,£ = 0), characterized by the zero-th order dis- 
persion uj = Vowo and by density oscillations of the form 
Sn ~ (r 2 — 3/5i? 2 ). In this case Eq.(|iT|) yields the rele- 
vant result 
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showing that the fractional shift of the monopole fre- 
quency is proportional to the square root of the gas pa- 
rameter evaluated in the center of the trap. This correc- 
tion exhibits the same dependence on the gas parameter 
as the quantum depletion of the condensate, although the 
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coefficient of proportionality slightly differs in the two 
cases. It is useful to write the gas parameter in terms of 
the relevant parameters of the system as B 



a 3 n(0) = 1^! ( N X/*J_ 

on \ a ho 
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where N is the number of atoms in the trap and aj w = 
(h/mivo) 1 ^ 2 is the oscillator length. Using, for example, 
N = 10 6 and a/ aho — 6 10 -3 , we predict a relative shift of 
1 %. A similar value is found for the quantum depletion 
of the condensate. Eq.([l3|) shows that in order to enhance 
the value of the gas parameter it is more effective to 
increase the value of the ratio a /aho rather than the value 
of N which enter the equation with a much lower power. 
In practice, however, it is not easy to obtain large values 
of a 3 n(0) and hence large frequency shifts. So far the 
achievement of high densities is in fact limited by three 
body recombinations which cause the atoms to leave the 
trap. 

The above shift of the monopole frequency should be 
compared with other corrections which might be rele- 
vant in actual experiments, like finite size, non-linearity 
and thermal effects. Finite size effects arise because even 
in the mean field scheme the Thomas-Fermi value y/Eu>o 
holds only in the large N limit. These corrections arise 
from the "quantum pressure term" in the equation for 
the velocity field which is ignored in Eq.(Q), and can 
be calculated by a proper perturbation procedure in the 
Gross-Pitaevskii equation. Using a sum rule approach 
one obtains the following result for the leading correc- 
tion to the monopole oscillation in the large N limit and 
isotropic trapping p3[ 



Sum _ 7 faho\ A -, f R 
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where C = 1.3 is a dimensionless parameter and R = 
ahoi^Na/aho) 1 ^ 5 is the radius of the system. For the 
surface quadrupole mode one finds that the fractional 
shift has opposite sign and is larger by a factor 5. 

It is worth noticing that the corrections to the Thomas- 
Fermi value due to non-mean ( fl2| ) and finite size (|l4|) ef- 
fects depend on different combinations of the relevant pa- 
rameters of the system. In fact finite size effects depend 
on the combination Na/a,ho, while non-mean field ef- 
fects depend on N 1 '^a/aho- I n the thermodynamic limit, 
where N — > oo and ujq — > 0, with the product Nu>q kept 
constant, finite size corrections go to zero, while the gas 
parameter (jl^) has a finite value (for a discussion of this 
limit in trapped Bose gases see, for example, Q). The 
large N, thermodynamic limit is reasonably well realized 
in experiments. For example, using the same values for 
N and a/aho employed above, one finds that the finite 
size shift (|lj) of the monopole frequency is much smaller 
(~ 0.1 percent) than the non-mean field correction (12). 



In general finite size effects are negligible if the condition 
N ^> (aho/ a) 2 log(R/aho) is satisfied. 

Non-linearity is another important effect to discuss. In 
fact in actual experiments the amplitude of the oscilla- 
tion cannot be made arbitrarily small. The effects of 
non-linearity have been investigated in details in H] in 
the framework of the Thomas-Fermi approximation. The 
leading corrections to the frequency shift can be written 
in the form 5ui/u> = A 2 5 where A is the fractional ampli- 
tude of the oscillation of the atomic cloud confined in the 
trap and the coefficient S can be calculated in an explicit 
way [O. For the monopole mode in the spherical trap 
one has S = —1/6 so that for fractional amplitudes less 
than 10 %, the effects of non-linearity are very small. 

In addition to finite size and non linearity effects, one 
should also take into account that experiments are car- 
ried out at finite temperature. At present there is no 
fully reliable theory to account for the temperature de- 
pendence of the collective frequencies of these trapped 
gases. However one expects that these effects should van- 
ish very rapidly when kT is smaller than the chemical 
potential. A rough estimate of the thermal effect can be 
obtained by assuming that the shift of the real part of 
uj is of the same order as its imaginary part which is re- 
sponsible for the damping of the oscillation. This might 
provide an experimental control of the thermal effect on 
the frequency shift also in the absence of accurate low 
temperature thermometry. 

Finally an important question concerns the role of the 
anisotropy of the confining potential. In fact most of 
magnetic trap are at present non spherical. Furthemore, 
high density values are at present more easily reached 
working with highly deformed traps. For an axially de- 
formed trap of the form V ext — \raio\r\ + \mJ\z 2 , 
(14) where rj_ = (x 2 + y 2 ) 1 ^ 2 is the radial coordinate, the 
HD equations (|9|) admit several interesting solutions. In 
addition to the dipole (center of mass) oscillation, the ex- 
citations so far investigated experimentally are the m = 2 
quadrupole mode, whose density varies as Sn ~ r 2 Yxi^ 
and the m = oscillations resulting from the coupling 
between the quadrupole and the monopole modes (notice 
that the z-th component, m, of angular momentum is still 
a good quantum number in axially deformed systems). 
The m = 2 quadrupole mode has frequency lu = 2uj±_ and 
is not afffected by non-mean field effects since V 2 <5n = 0. 
Conversely the decoupled frequencies of the m = modes 
are given by Q 
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where A = u z /lu± characterizes the asymmetry of the 
trap. The corresponding density oscillations have the 
form 5n — (-2(i(s 2 - 2)/muj 2 + r\ + (s 2 - A)z 2 ). After 
some length, but straightforward algebra, one derives the 
result 
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for the frequency shift of the m = modes where 



(16) 



/±(A) = I 
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and the index ± refers to the higher (+) and lower (— ) 
solutions of Eq. (|l5|) . Notice that for a spherical trap 
(A = 1) the solutions of (Jl^) are s 2 = 5 (monopole) and 
s 2 = 2 (quadrupole) . In the first case one finds f + = 1 
and one recovers result (|I^) for the breathing mode, while 
in the second case there is no shift. In the figure we show 
the functions /+ and /_ relative to the two modes as a 
function of the deformation parameter A. 

Another important consequence of the use of deformed 
trap concerns the effects of non-linearity. In fact it has 
been shown G3| that these effects can be amplified or 
reduced by changing the value of A. For example, choos- 
ing A = 1.40 one finds that the non-linearity coefficient 
S of the high lying solution of ( |l5| ) vanishes. For this 
mode one finds s 2 = 7.1 and / + = 0.88. Another inter- 
esting case is obtained choosing A = In this case the 
non-linear effect is vanishingly small for the low energy 
solution. For this mode one finds s 2 = 3.2 and /_ = 0.38. 

As already pointed out the largest values of the gas 
parameter have been so far reached for cigar trap con- 
figurations. In this case the frequency shift of the lower 
mode with dispersion s 2 = 5A 2 /2 is reduced significantly 
(/_(0) = 1/6) with respect to the case of the A = 1 
breathing mode (/+(1) = 1). Conversely the higher 
mode, which corresponds to a compressional radial os- 
cillation with dispersion s 2 = 4, exhibits only a small 
reduction (/+(0) = 5/6). It is also worth pointing out 
that for this mode the non-linear effect vanishes as Jli| ] 
8lo/uo — — 0.0938A 2 A 2 . Finite size effects are also van- 
ishing in the same limit. The fact that non linear and 
finite size effects vanish for A — > reflects the occurrence 
of a hidden symmetry [19] of the Gross-Pitaevskii equa- 
tion characterizing the radial compressional mode in two- 
dimensional gases and migth be used to improve the ac- 
curacy of experimental measurements. One should how- 
ever notice that, for very small values of A, the trapped 
gases exhibit, in addition to the radial excitations, also 
a low- lying branch of axial modes p6[ . This could pro- 
duce a "parametric" instability jlTf of the compressional 
radial oscillation due to decay into two or more axial ex- 
citations. 

In conclusion in this letter we have derived the first 
corrections to the collective frequencies of trapped Bosc 
gases arising from effects beyond mean field theory. We 
have shown that with reasonable choices of the parame- 
ters these effects, although small, might be visible exper- 
imenatally. Their direct observation would represent an 
important achievement in the study of many-body effects 
associated with the occurrence of Bose-Einstein conden- 
staion. 
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Figure Caption. Functions /+ and /_ relative to the 
higher and lower m = modes (see eqs.(|i"5|) and ([To])), 
as a function of the deformation parameter A = lu z /ujj_. 
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